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Abstract 

We investigate baryogenesis in the z^MSM (neutrino Minimal Standard Model), which 
is the MSM extended by three right-handed neutrinos with masses below the electroweak 
scale. The baryon asymmetry of the universe can be generated by the mechanism via flavor 
oscillation of right-handed (sterile) neutrinos which are responsible to masses of active 
neutrinos confirmed by various experiments. We present the kinetic equations for the 
matrix of densities of leptons which describe the generation of asymmetries. Especially, the 
momentum dependence of the matrix of densities is taken into account. By solving these 
equations numerically, it is found that the momentum distribution is significantly distorted 
from the equilibrium one, since the production for the modes with lower momenta i<CT 
(T is the temperature of the universe) is enhanced, while suppressed for higher modes. As 
a result, the most important mode for the yields of sterile neutrinos as well as the baryon 
asymmetry is k ~ 2T, which is smaller than (k) inferred from the thermal average. The 
comparison with the previous works is also discussed. 



1 Introduction 



The baryon asymmetry of the universe (BAU) is one of the most puzzling problems in particle 
physics and cosmology. It is certain that the Minimal Standard Model (MSM) cannot account 
for its origin. There have so far been proposed various scenarios of baryogenesis by considering 
physics beyond the MSM. (See, for example, a recent review pp.) Among them leptogenesis j2] 
by superheavy right-handed neutrinos is one of the most motivated scenarios. This is because 
these new particles allow us to give non-zero masses to neutrinos which have been confirmed 
in many oscillation experiments. Being singlet under the MSM gauge group they can obtain 
Majorana masses which are completely independent on the electroweak scale. Right-handed 
neutrinos having superheavy masses then can generate the lepton asymmetry by their decays 
which can be a source of the BAU. In addition, the observed smallness of active neutrino masses 
is naturally explained by such fermions through the seesaw mechanism [3] . The required masses 
are so heavy that it is impossible to directly test these new particles in near future experiments. 

Akhmedov, Rubakov and Smirnov have proposed j3] another attractive scenario of baryo- 
genesis by using right-handed neutrinos (mechanism via neutrino oscillation). One of the most 
important features is that the mechanism works when their masses are smaller than the elec- 
troweak scale, which are within the reach of current experiments. The flavor oscillation among 
right-handed neutrinos in the early universe leads to the separation of the lepton number into 
right-handed neutrinos and left-handed leptons. The sphaleron process then converts the asym- 
metry of the left-handed leptons into the baryon asymmetry for high temperatures [5]. 

This mechanism can be realized in a simple and attractive framework called as the neutrino 
MSM (z/MSM) [6j [7J in which three right-handed neutrinos are introduced with Majorana 
masses below the electroweak scale. Because of the very suppressed neutrino Yukawa interaction 
the seesaw mechanism is still effective, which leads to three active neutrinos Vi (i — 1, 2, 3) and 
three sterile neutrinos Nj (I = 1,2,3) as mass eigenstates. The former ones are responsible 
to the neutrino oscillations observed in experiments. The latter ones give the solutions to the 
cosmological problems. The lightest sterile neutrino Ni with mass ~ 10 keV can be a candidate 
of dark matter, while the rest two N2 and with masses ~ 1 GeV can generate the BAU 
through the mechanism above. (For a review see Ref. [8].) Furthermore, the non-minimal 
coupling of Higgs field to gravity allows to realize the cosmic inflation [9] . 

Various aspects of baryogenesis in the z/MSM has been studied until now [7J QUI dH IT2"] . 
The generation of the asymmetry is described by the matrix of densities [131 [EH H5] of sterile 
(right-handed) neutrinos and active (left-handed) leptons. The previous works were based on 
the kinetic equations for the matrix of densities in Ref. [TJ, where the terms describing the 
transfer of asymmetries between sterile neutrinos and active leptons are added to the original 
ones in [4j. As shown in [TJ, such terms are essential to generate enough amount of the BAU 
in the z/MSM. In spite of their significance such terms are introduced by a heuristic approach. 

The main motivation of this article is to improve the estimation of the BAU in the z/MSM. 
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One of the most unsatisfactory points in the previous works is that the evolution of the matrix of 
densities has been traced only by an approximate way. Namely, the kinetic equations analyzed 
previously are for the typical, single mode with momentum k ~ T (T is the temperature 
of the universe), which is expected to give a dominant contribution to the asymmetry. By 
extrapolating the result of such a mode to other modes the asymmetry in the number density 
is estimated. 

In this article, thus, we derive the kinetic equations in which the momentum dependence 
of the matrix of densities is fully taken into account. For this purpose we shall re-evaluate the 
destruction and production rates of sterile neutrinos as well as active leptons paying attention 
to the two points. One is the momentum dependence of the particle of interest. The other is 
the corrections coming from the fact that the states in the destruction and production processes 
may differ from the thermal equilibrium states (i. e., the deviations from the thermal equilibrium 
for sterile neutrinos and the chemical potentials for active leptons). As we will show later, the 
kinetic equations are written as simultaneous integrodifferential equations for the matrices of 
densities of sterile neutrinos and active leptons. 

Interestingly, our kinetic equations include the terms connecting between sterile and active 
sectors, which is crucial for the baryogenesis in the z/MSM as mentioned before. It will be 
shown that such terms arise automatically as the corrections to the destruction and production 
rates owing to the deviations from the equilibrium states in the scattering processes. It turns 
out that the coefficients of these terms are different from [7] and, in addition, there appears 
a new type of terms which couples sterile neutrinos to active anti-leptons (rather than active 
leptons). 

This paper is organized as follows: In Sec. |5]we briefly review the framework of the analysis 
and the baryogenesis scenario in the z/MSM. In the Sec. |3j we derive the kinetic equations 
which describe the generation of the baryon asymmetry paying attention to the momentum 
dependence in the matrices of densities. Sec. @] is devoted to study the numerical solutions to 
the obtained equations. We show the mode by mode evolution of the matrices of densities and 
their momentum distributions. The comparison with the previous works is done in Sec. We 
perform the quantitative comparison of the yields of the BAU, but also clarify the differences 
between the kinetic equations in the literature. We conclude in Sec. El 

2 Baryogenesis in the z/MSM 

At the beginning we review the z/MSM, which is the MSM extended by three right-handed 
neutrinos v^i (I = 1, 2, 3) with Lagrangian 

£ z/Msm = £ M sm + iJ7 Ril fl d» v Ri ~ (Fad L a $ vri + -y- v R \ v RI + h.c.^j . (1) 

where £msm is the MSM Lagrangian, F aI are neutrino Yukawa coupling constants, and $ 
and L a (a = e, u, r) are Higgs and lepton weak- doublets, respectively. The Majorana masses of 



right-handed neutrinos are denoted by Mj. After the electroweak symmetry breaking neutrinos 
also obtain the Dirac masses, [Mjj] a ; = (($) is a vacuum expectation value of the Higgs 

field), and tiny neutrino masses can be explained by the seesaw mechanism if |[Md] q j| <C Mj is 
satisfied. In this case, mass eigenstates of neutrinos are consist of active neutrinos z/j (i — 1, 2, 3) 
and sterile neutrinos Ni (I = 1,2,3). Active neutrinos are the mass eigenstates of ordinary 
neutrinos and sterile neutrinos are almost the right-handed states Ni ~ Uri with masses given 
by the Major ana masses Mj approximately. 

On the other hand, when we discuss baryogenesis, the temperatures of interest is higher 
than the electroweak scale, and we can neglect all the masses of the MSM particles and also 
the mixing between active and sterile neutrinos. In this case we can treat active and sterile 
leptons as left-handed leptons and right-handed neutrinos, respectively. 

Among three sterile neutrinos, the lightest one N\ is a candidate of dark matter. In this dark 
matter scenario, the coupling constants of neutrino Yukawa interaction of N\ are required to be 
very small. (See, for example, a review [§].) As a result, N\ gives the negligible contributions to 
the seesaw mass matrix of active neutrinos as well as the BAU, and then two sterile neutrinos 
N 2 and N 3 are responsible to these two phenomena. 

In the z/MSM the sufficient amount of baryon asymmetry can be generated [7J by the 
mechanism via neutrino oscillation [I]. The flavor oscillation between sterile (right-handed) 
neutrinos N 2 and N3 are induced by the medium effect, which can generate the asymmetries 
in leptons together with the CP violation in neutrino Yukawa matrix. The asymmetry in 
active (left-handed) leptons is then transformed into the baryon asymmetry due to the B + L 
breaking sphaleron process for high temperatures T > T w = C(10 2 ) GeV [5]. Notice that 
sterile neutrinos is possible to be out of thermal equilibrium due to the smallness of neutrino 
Yukawa coupling constants. If this is the case, all the conditions for successful baryogenesis [16] 
can be satisfied. 

In order to describe a series of the processes generating asymmetries, we have to deal with 
the coherent evolution of sterile neutrinos with the flavor oscillation, and also the incoherent 
scattering processes with surrounding medium for the destruction and production of sterile 
neutrinos. We then use the formulation with the matrix of densities [121 [HI [32] . The matrix of 
interest here is one in the flavor space consisting of sterile neutrinos N 2 and N 3 with positive 
helicities, N 2 and N 3 with negative helicities, active leptons L e , and L T , and their anti- 
particles L e , and L T , which leads to a 10 x 10 matrix. The diagonal elements of this 
matrix are nothing but the usual occupation numbers, and the off-diagonal elements contain 
correlations of the flavor mixings. 

In the situation under consideration, however, it can be simplified as described in Ref. [7J. 
First of all, the conservation of the total lepton number holds for the temperatures of baryo- 
genesis (T > TV), since Majorana masses Mj for sterile neutrinos are sufficiently small in the 
z/MSM. We can thus neglect all the elements which break the lepton number. Second, active 
leptons possess gauge interactions and Yukawa interactions with right-handed charged leptons 
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in addition to neutrino Yukawa interactions, which induce large energy gaps between active 
and sterile states through the thermal effect. Then, the transitions between active and sterile 
leptons are highly suppressed, and the elements corresponding to the mixing between L a and 
Ni (and also L a and Nj) are neglected. At this stage, the matrix of densities in the system is de- 
composed into two 2x2 matrices pn and p^ for and iV^, and two 3x3 matrices pi and pi 
for active leptons L e ^ T and L e ^ T . As in Ref. [TT] . pi is considered to be pi = p u + p e = N D p u , 
i.e., the sum of two contributions of lepton doublet. Finally, the off-diagonal elements of p^ 
and pi can be neglected for the temperatures of interest. This is because the flavor transitions 
among active leptons are also suppressed due to the medium effects induced by the Yukawa in- 
teractions of charged leptons with hierarchical coupling constants. Note that active leptons can 
maintain the kinetic equilibrium due to the rapid interactions with medium, and hence their 
diagonal matrices of densities can be expressed by using the dimensionless chemical potentials 
p Ua (the ordinary chemical potential divided by temperature) as 

PL (k) = N Df f*(k) A , PL {k) = N D p eq (k) A- 1 , (2) 

where No = 2 is an SU(2) factor, A = diag(e M " e , e^"' 1 , e M " r ) and p eq (k) is the equilibrium 
distribution function. Throughout this analysis we apply the approximation of the Boltzmann 
statistics in which p eq (k) = e~ k l T . 

The evolution of the system is now described by two 2x2 matrices of densities pN{k) and 
Pft(k), and three differences of the distribution functions [piik) — pi(k)] aa (or three chemical 
potentials fi Va ). The previous works on baryogenesis via neutrino oscillations have been based 
on the kinetic equations for these variables given in Ref. [1] or [7j. The original work [3] has 
considered the kinetic equations for p^ and p^ without including the effects from active leptons. 
As pointed out in Ref. [TJ, their equations give too small baryon asymmetry to account for the 
observed value in the t/MSM. This is because the lightest sterile neutrino N± is a dark matter 
candidate and essentially only two sterile neutrinos N 2 and A^ 3 participate in baryogenesis. 
Furthermore, Ref. [7j has introduced the terms which connect active and sterile sectors, and 
solved the kinetic equations for both p N $ and p L g. It has been shown that the additional terms 
boost the generation of the asymmetry and hence are vital for the successful baryogenesis in 
the z/MSM. 

There are, however, various aspects to be improved in the previous works towards the 
precise prediction of the baryon asymmetry in the i/MSM. One of the most important points is 
to include correctly the momentum dependence in the matrices of densities pN{k) and pfj{k). So 
far, it has been analyzed only the evolution of the single mode with a typical momentum k ~ T, 
and estimated the asymmetries in number densities by the approximation that the occupation 
numbers are proportional to the equilibrium one p eq (k). To get rid of this uncertainty, therefore, 
we will write down the kinetic equations of the matrices of densities for all the modes and solve 
them numerically taking care of the expansion of the universe. It will be shown later that the 
momentum distributions of p N ^ are distorted significantly from p cq . 
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Figure 1: Feynman diagram for self energy of sterile neutrinos. 

3 Kinetic Equations 

Now we are at the position to derive the kinetic equations for p N ^ and p L i. First of all, let 
us consider the time evolution of p^. Our construction is based on Ref. [15] (as in (H [7]) and 
starts with 

= -i[H N (k N ), PN (k N )} - \{r d N (k N ), PN (k N )} + ^{r p N (k N ), i - PN (k N )} , (3) 



where 1 denotes the unit matrix with appropriate dimensions!^!] Hn is the effective Hamilto- 
nian, H N = H N + Vjv, where the free part is [H N {k N )]u = E Ni 5jj with E Nl = \Jk 2 N + Mj and 
Vjv is the effective potential induced by the medium effects. Y d N and Y P N are the destruction 
and production rates of Nj. From now on we shall apply the approximation of the Boltzmann 
statistics and replace the third term of Eq. (j3]) as \{Y p Nl 1 — Pn} —> T^. 

The first term of Eq. ([3]) describes the coherent evolution of p^ and the oscillation of sterile 
neutrinos occurs due to the off-diagonal elements of Vjv, which is essential for baryogenesis 
under consideration. It is found from the self energy for sterile neutrinos at finite temperatures 
in Fig. [T]that the effective potential for the mode k = k^ is given by [T9] 

N T 2 



#2 



where we disregard the correction to Vn from the asymmetries in active leptons, 

In the estimation of Vn (as well as r^ p below) all masses including Mj are neglected since 
they are irrelevant for temperatures of interest. (Note, however, that we keep Mi in H N because 
they are crucial for the oscillation of sterile neutrinos.) Further, we first calculate them in the 
basis where neutrino Yukawa matrix is diagonal, and then find the expression in the original 
basis shown in Eq. ([1]). 

Let us then estimate the destruction and production rates of Ni with momentum k^- In 
the considering temperatures the dominant contributions come from the scattering processes 
(A) iVj + Q L <-> L a + t R , (B) Ni + i R <-> L a + Q L , and (C) N T + L a ^ t R + Q L [1], shown in 



* x We have neglected the non- linear effects of pm since the interaction rates between sterile neutrinos are 
sufficiently small. Otherwise, see Ref. [17l H8] . 

# 2 We have numerically confirmed that the change of the final baryon asymmetry by this effect is negligibly 
small. 
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Figure 2: Feynman diagrams for scattering processes of production and destruction rates. 

Fig. [2J Here Ql and tn denote left-handed quark doublet of third generation and right-handed 
top quark. We then divide the rates into three parts: 



r d N p (k N ) = rf w (k N ) + rf( B \k N ) + rjf (c) (^) 

The destruction rates of each process are found to be 

[r# A) (*")L = [rTikN)}^ = j d N (k N ) [F*F] 



(5) 



i j 



[T d f \k N )] IJ = 7 d N (k N ) [F^F] IJ +[6Tj f (k N )] 
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Here we have introduced 



7&(* 



N c N D h 2 t T 2 



64vr 3 k N ' 

where Nq = 3 is a color factor and ht — 1 is the top Yukawa coupling constant, and 



(6) 



(7) 



[5T d N (k N )] IJ = Y N (k N ) 



dk L k L 



F\p T L {k L )-N D p c \k L )i)F 



J i J 



/o N D T 2 

= 1 d N (k N )[F^A- 1 -l)F] IJ} (8; 
where we have used Eq. ([2]) in the last equality. On the other hand, the production rates are 
K {A \kN)] u = [rT\k N )] u = 7 d N (k N )r(k N ) [F^F] U + [6T p N (k N 



i j ' 



K (c \k N )] „ = i d N (k N ) r(k N ) [f*f] u , 



(9) 



where 



f°° dk L 1 - P ^(k N ) 



l kN N D k N p^(k N ] 
■■l d N {k N ) P c \k N ) [F\A-1)F]^ 



F\p L {k L )-N D p^{k L )l)F 
F^(p L (k L )-N D p^(k L )l)F 



u 



(10) 



It should be noted that STf^ and 5T P N vanish when all the chemical potentials of active 
leptons become zero (i.e., piik) = pi(k) = Ncp eq (k)l). In this case, the rates of each process 
satisfy I^y A ' B (&jv) = p QCL (k N )r d N A ' B ' C \kN) , and the total rates are given by 

K eq (^)] 7J = 3<y£(**r) [^F] u , [T^(k N )] = p cq (k N ) [T d N e(l (k N )] IJ . (11) 

The kinetic equations for p^ and p^ are then summarized as 



df>N ^ N) = -i [H° N (k N ) + V N (k N ), p N (k N )} - l - {T d N c \k N ), p N (k N ) - p c \k N )\] 

+25T p N (k N ) - l - {5T d N (k N ), p N (k N )} . (12) 

The equation for p^ can be found by the CP conjugation, i.e., by exchanging F -H- F* and 
Pl PL (or A A' 1 ). 

Next, we turn to consider the left-handed leptons. Similar to the above case, we start with 
the equation for p L as 

= -i [H u (k L ), PL (k L )] - \ {T d u (k L ), p L (k L )} + N D Tl(k L ) . (13) 

Here remember that pi denotes the sum of SU(2) doublet contributions pi = p u + p e . As 
explained in Eq. fl2]), only the diagonal elements of both sides are taken into account and then 
the first term of right-hand side becomes irrelevant for the discussion. Further, active leptons 
maintain the kinetic equilibrium as Eq. fl2]) due to the large interaction rates, we have only to 
consider the evolution of their chemical potentials. 

«- h ^ = ik/^sW fe) -« (fe) ]-- (14) 

We then estimate the destruction and production rates of active (left-handed) neutrinos, Tf, 
and r£, caused by neutrino Yukawa interaction. Similar to sterile neutrinos, the dominant 
contributions are found from the scattering processes (A), (B) and (C). The production rates 
are then estimated as 

[r; (C, fe)]„, =lt(k L )[FF\ fl + [5Tt(k L )} afr (15) 



where 



= (16) 
-^[F{pl{k N )-r{k N )l)F^] ap . (17) 

On the other hand, the production rates for three processes are 

[n {A \k L )] a p = [T p J B \k L )]^ = H{h L ) p e \k L ) [FF\ p + [5Vl{k L )] ap , 
[K iC) (k L )} a p = li{k L )r{k L ) [FF*]^, (18) 
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where 



WM^ = HIV P m (^ \ /* TT 1 -^ 1 [*W*»> - „«-(Mi)f ' 




a/3 



fcL &z P^l) 



F(p 7V (^)-p eq (A: J v)l)Ft L (19) 

J a/3 



The rates of left-handed charged leptons are the same as those of active neutrinos, and the 
CP conjugation (F -h- F* and p^ 4-> p^r) gives the rates for v. The equations for chemical 
potentials are then given by 



d J^ = - 1 i{T) [FFtj^tanh^ 

- tanhp^ [Fp JV (A; JV )Ft + F*p R (k N )F T ] a Jj . (20) 

Therefore, the three equations, f lT2"j) for p^v, the CP conjugation of (fT2"j) for p^, and (J2"0]) for 
p Va , are the kinetic equations for our study of baryogenesis in the i/MSM. 

Here are some comments. First, we have fully taken into account the momentum depen- 
dence in the matrices of densities and p^. As a result, the kinetic equations are written as 
simultaneous integrodifferential equations, which allow us to find the distribution of the occu- 
pation numbers of sterile neutrinos, the significant modes to generate the baryon asymmetry, 
and the temperatures crucial for baryogenesis. 

Second, we have calculated the destruction and production rates of sterile neutrinos and 
active leptons for a given momentum. It should be noted that we have taken into account 
the case when the initial and final states in scattering processes which are not in the thermal 
equilibrium and estimate the corrections to the rates from the deviations p N jq — p eq and p L i — 
Nr>p eq - Interestingly, the communication terms between p^ and Pl (pr an d Pi) introduced in 
Ref. [7] arise automatically, which is the term with ST P N in Eq. ( fl2|) . Now we can understand 
well the origin of such terms. For instance, 5T P N originates in the correction to the production 
rate of iV in the scattering processes (A) and (B) when chemical potentials of active leptons 
are non-zero. The similar discussion can be applied to the communication terms in Eq. (120 p . 
As shown in Ref. [7], they play a crucial role to generate the BAU in the z/MSM. 

Finally, our kinetic equations also contain the communication terms between p^ and pi (p^r 
and pi) which have not been discussed before. For example, the term with 5T d N in Eq. (1121) 
corresponds to it. Such a term arises as the corrections to the destruction rates of sterile 
neutrinos due to the non-zero chemical potentials of active leptons via the scattering process 
(C). The important point is that the term with 5Y d N is the second order of the matrices of 
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densities (i.e., PnPl) an d the kinetic equations are no longer linear in Pn,Pn an d fJ> Va - Further, 
the presence of these terms is essential to ensure the conservation of the lepton number. 

The total lepton number becomes a conserved charge for temperatures of baryogenesis since 
the Majorana masses are much smaller than Tw, as mentioned before. In the considering 
system, this leads to 



where R is the scale factor of the expanding universe and nx is the number density of particle 
X. It can be written in terms of the matrices of densities (and chemical potentials) as 



A simple calculation shows that our equations (TT2T) and (1201) satisfy this equality. It is found 
for the scattering processes (A) and (B) that the production (or destruction) terms in p N 
and cancel with the destruction (or production) terms in p L and pi, respectively, after 
taking the trace and integrating over momentum. As for the process (C) the production (or 
destruction) terms in p^ and p^ cancel with the production (or destruction) terms in pi and 
Pl, respectively. Thus, if the system starts with the lepton symmetric universe, the asymmetry 
of sterile neutrinos is always opposite to that of active leptons. Namely, the kinetic equations 
here describe not the generation of the total lepton asymmetry (i.e., leptogenesis), but the 
separation into sterile and active sectors. 

4 Numerical Solution of Kinetic Equations 

In this section, we shall study the numerical solution of the kinetic equations ffT21 and (|2"U|) . 
To incorporate the expansion of the universe we replace the time derivative in Eq. (|T2|) as 
f t g - Hk-j^, where H = T 2 /M (M = 7.12 x 10 17 GeV) is the Hubble expansion rate. 
The time-temperature relation is then given by ^ = j^. In this analysis, we take the initial 
conditions as pN(k) = Pfi(k) = p Ua = 0, which may be realized by the primordial inflation. We 
solve these equations till T = T w in order to estimate the BAU, where T w is the sphaleron 
freezing temperature. (See, for example Ref. |2T].) Following to Ref. [12] , we set Tw = 140 
GeV from now on. 

In solving the kinetic equations, we have to specify neutrino Yukawa coupling constants 
F a j and masses of sterile neutrinos M 2 and M 3 . The estimation of the baryon asymmetry 
in full parameter region is beyond the scope of this analysis, which is postponed in another 
publication. Here we choose a specific choice of parameters and study the properties of the 
solutions. 
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Figure 3: Evolution of [pN(k)] 2 2 with k/T = 0.02, 0.1, 0.33, 1, 2 and 4-4 ^ n the left panel, and 
evolution of Y/v 2 3 in the right panel. 



As for the neutrino Yukawa couplings we follow the notation in Ref. [22]. In this analysis 
we consider the normal hierarchy of active neutrino masses and take m 3 = 4.89 x 10~ 2 eV, 
m 2 = 8.71 x KT 3 eV, m x = 0, sin 2 6 23 = 0.42, sin 2 6 12 = 0.312 and sin 2 13 = 0.025 [20]. Further, 
we choose the parameters of sterile neutrinos as£ = +l, w = 7r/4, 5 = 7n/4: and rj = n/3. We 
write masses of sterile neutrinos as M 3 = Mjv + AM^/2 and M 2 = — AM^/2, and choose, 
as representative values, = 10 GeV by fixing AM^/M^ = 10~ 8 otherwise stated. 

The diagonal elements of the matrices of densities pn and p?q are the occupation numbers 
of Nj and Nj. The evolution of these quantities are shown in the left panel of Fig. [3j In the 
present choice of parameters (i.e., u = 7r/4) the occupation numbers of N 2 and N 3 are almost 
the same, and we only show the evolution of [ptv]22- It is seen that the occupation numbers 
increase as the temperature decreases due to the production by the scattering processes (A), 
(B) and (C). In Eq. (fT2"|) the third term on the right-hand side which proportional to p eq is 
the main source of the production of sterile neutrinos. Since the destruction rate in Eq. fill I) 
is inversely proportional to the momentum, the production of modes with lower momenta is 
more effective. We find that the low modes stop to grow eventually and take constant values 
given by [pN(k)]n = [pN(k)]n = p eq (k) afterward. This thermalization temperature becomes 
higher for the mode with smaller momentum. 

The yields of Nj and Nj, defined by the ratio between the number and entropy densities, 
are given by 

^ [ P N(k)] H , Y r 1 r d " k 



Y, 



(23) 



(2tt) 3 Lnn /J «' JV/ sj (2tt) 3 

where the entropy density is s = ^-g s T 3 with g s = 106.75. When the system is fully thermalized 
(i.e., [pN(k)} = [pffik)] = p eq (k)l), the yields take a constant value 
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Figure 4: Evolution o/Re[pjv]23 (fed solid lines) and Im[pjv]23 (blue dashed lines) for the modes 
with k/T = 0.1 (left panel), 1 (middle panel) and 10 (right panel), respectively. 



for the temperatures of interest. It is found that Y^ 2 and Y^ 3 increase as 1/T until they close 
to the equilibrium values as shown in Fig. |3j 

On the other hand, the off-diagonal elements of the matrices of densities contain correlations 
of the flavor mixing. In Fig. H]we show the evolution of 23 for the modes with k/T = 0.1, 

1 and 10. It is seen that both real and imaginary parts of [pN{k)]23 start to oscillate around 
the temperature [7] 



T osc (k) 



M A.U \ -\/ 



A? 



3 (k/T) 



1/3 



(25) 



and hence the oscillation for the mode with lower momentum begins at higher temperature. 
We should note that the amplitude of the oscillation becomes damped and the off-diagonal 
elements of pN vanish when the mode gets in the thermal equilibrium. This behavior can be 
seen for the mode with k/T = 0.1 in Fig. HI 

The flavor oscillation between N2 and N3 can generate the asymmetries of sterile neutrinos 
as well as active leptons [HH], together with the CP violation in neutrino sector. The analytical 
description of these production processes are found in Refs. jU [TJ QUI E], and we just show the 
results of the numerical analysis. 

The asymmetries in the occupation numbers of sterile neutrinos, [Ap N (k)}n = [pN(k)]n — 
[pN{k)}ni are generated as shown in the left panel of Fig. [5j It is seen that such an asymmetry 
for the mode with k is generated at T ~ T osc (k) and the absolute value of the asymmetry 
becomes larger for the lower modes. We also find that the asymmetries of the mode with 
k/T = 0.02 begins to damp at Tw/T ~ 0.2, at which the mode gets in thermal equilibrium. 
Having the asymmetries of each mode, we can estimate the asymmetry in the number density 
as 
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(26) 
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Figure 5: Evolution of [Ap^]22 (solid lines) and [Apn]33 (dotted lines) with momenta k/T = 
0.02, 0.1, 0.33 and 1 in the left panel, and evolution of Y&n 2 , Yajv 3 an d ^AiVtot ^ n ^ e r ^9^ 
panel. 
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Figure 6: Evolution ofY^ e , Yal andY^L T in the left panel, and evolution ofY&L tot andY^Nu 
in the right panel. 



As shown in the right panel of Fig. \5\ Faa^ increases as ~ 1/T as long as the system is away 
from the equilibrium state. We also show the evolution of total asymmetry of sterile neutrinos, 

^ATVtot = Yan 2 + ^AAT 3 ■ 

In the considering system, the non-zero asymmetry in the sterile sector necessarily indicates 
the presence of asymmetry in active sector. The yields of active flavor asymmetries are expressed 
in terms of chemical potentials as 

Y/\L a = t~ smh p a . (27) 

The evolution of Y^L a is shown in the left panel of Fig. EJ The asymmetries of active flavors Yal 
are larger than those of sterile neutrinos Faatj- This is because Yal q is induced at the fourth 
order of neutrino Yukawa coupling constants, while Y& Nl is at the sixth order [7]. As shown in 
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Ref. [IT] when Imu; = 0, the yield Val q is proportional to the CP asymmetry parameter A^ 2 . 
In the present choice of parameters, they are 

^32 : ^32^32 = -1-0: -15: 16. (28) 

It is then found that our numerical solutions of the kinetic equations reproduce the analytic 
relations of asymmetries of active lepton flavors. 

The total asymmetry stored in active leptons is now obtained as ^ALtot = + ^al„ + 
Yal t , which evolution is shown in the right panel of Fig. [6J For comparison we also show the 
asymmetry in sterile sector. It is clearly seen that Y^ tot + Y& Ltot = due to the conservation 
of the total lepton number in the considering system. These total asymmetries are obtained as 
the 0(F 6 ) effects. We also find that the generation of the total asymmetries becomes effective 
at the temperature 

T L ~T osc (k = T). (29) 

We find Ti ~ 6.2 x 10 3 GeV for our choice of parameters. These yields scale as 1/T, which are 
consistent with the estimation in Ref. [7J, even if one includes the momentum dependence of 
the matrices of densities properly. 

The asymmetry in active sector can be partially converted into the baryon asymmetry as 
AB = — ||AL tot due to the rapid sphaleron transition [5]. The BAU observed in the present 
universe is then given by 

Y B = ~ Y ALt J T=Tw = -1.53 x 10- 3 N D J2 sinh Ml J T=Tw . (30) 

a 

Notice that the asymmetry for T < T\y cannot contribute to the BAU since the sphaleron 
process is ineffective for the lower temperatures. In the considering parameter choice, the BAU 
is given by (See also the case III in Tab. [TJ) 

2.9 x 10- 10 for M N = 10 GeV 
*u=< , . ^ r . f " * (31) 



5.4 x 10" 14 for M N = 100 MeV 



When the sphaleron process is frozen at T = T\y, the distribution of the occupation number 
[Pn (k)]22 — [pN(k)} 33 is represented in terms of momentum in the left panel of Fig. [7J We should 
stress here that the momentum dependence in the occupation number [pn{^)]ii is significantly 
different from the equilibrium one p eq (k). This is because of the momentum dependence in the 
destruction and production rates of sterile neutrinos estimated in the previous section. These 
rates are characterized by ^%{k) in Eq. (JTj) and 'jff(k) oc 1/k. This results in the enhancement 
and the suppression of the production of the mode with lower and higher momentum than about 
T, respectively. Especially, the low modes with k/T < 0.1 are equilibrated until T = T w and 
coincides with p eq (k). It is also found that the most significant mode for the number density is 
found to be k ~ 2T, which is smaller than (k) = 3T obtained by the thermal average. 




Figure 7: Momentum distributions of [pN{k)}22 (red solid line) and p eq (k) (black dotted line) 
in the left panel. Momentum distributions of [Ap^(k)}22 (green dashed line), [Apjv(fc)] 3 3 (blue 
dotted line), and [Ap N (k)] tot (red solid line) in the right panel. The distributions are evaluated 
at T = Tw and are multiplied by (k/T) 3 . 

In the right panel of Fig. [7] we show the momentum distributions of the asymmetries of the 
occupation numbers, [Apjv(£;)] 2 2, [Ap^^k)^ and their sum [A/?7v(&)]tot> at T = Ty/- Similar to 
the occupation numbers, they have different momentum distributions from p eq (k). (See Fig. [HJ) 
Interestingly, the most important mode for the total asymmetry in sterile sector (and hence for 
the BAU) is found to be k ~ 2T as before. 

Therefore, these characteristic behaviors of the distributions show that the momentum 
dependence in the matrices of densities (as well as the destruction and production rates) is 
indispensable for the correct evaluation of the BAU in the z/MSM. 

5 Comparison with previous works 

In this section we shall compare the results obtained in this analysis with the previous works, 
especially given in Ref. [11], and also estimate the qualitative difference on the BAU. So far the 
momentum dependence of p^ and p^ have been approximately taken into account. To make 
the differences clear, let us first consider what kind of approximations should be imposed on 
the kinetic equations (TT2"j) and (|2"U|) in order to reproduce the previous results. 

It is then found that there are three important assumptions, (i) The momentum dependence 
is approximated by the equilibrium one, namely, 

[p N (k)] ZJ = [R N ] IjP ^(k) , [ P N(k)] u = [Rx] IjP e \k) , (32) 

and the evolution of Rn and R^ is studied, (ii) The thermal average is taken for the destruction 
rates of sterile neutrinos, (iii) The destruction and production rates by the scattering process 
(C) are identified with those by the process (A) or (B). We will discuss these points in turn. 



1 A 



First of all, based on the assumption (i), the equations for and chemical potentials are 
from Eqs. ^ and ([20]) as 



dR 



N 



dt 



dt 



-i[H N (k), R N ] 



{F*F, R N - 1} +2^) F\A - 1)F 



2 

37, d (T) 



{F^A -1 - 1)F, R N ] 



(33) 



[FFt] aQ tanh// a + ?£21[FR N F* - F*R R F T ] i 



' cosh fj, c 



{ [F(R N - l)Ft] QQ (l - tanh^) - [F*(R N - l)F T ] aQ (l + tanh^)} .(34) 



The equation for R^ can be found as explained in Sec. [3j We then obtain the coupled differential 
equations rather than the coupled integrodifferential equations. Let us recall here the origins 
of the communication terms in these equations. As for Eq. (1331) the third term comes from 
5T P N by the processes (A) and (B) and the fourth term comes from STfj by the process (C). 
Similarly, 5T^ by (A) and (B) induces the third term and STf by (C) gives the fourth term in 
Eq. (El. 

We can see that, although R N is introduced to be independent on the momentum k, the 
equation fl33l does depend on k through r yf [ (k) to which all the destruction and production 
rates of sterile neutrinos are proportional, and also through the effective Hamiltonian H^(k). 
Note that the rates of chemical potentials are independent on the momentum since they are 
proportional to 7^(T). In order to avoid this difficulty one might consider that the momentum 
in Eq. (1331) is replaced as k = k* ~ T, since such a mode is expected to give the dominant 
contribution to the number densities and also the asymmetries. However, we have to carefully 
choose k* to keep the conservation of the total lepton number. Under the considering situation, 
Eq. fl22|) is satisfied when 



= tr 



dR]\[ 




dRff 




dt 


fc=fc* 


dt 


k=k, 



d±_ N dA- 1 
D dt D dt 



(35) 



It can be shown that it is fulfilled if and only if k* = 2T. 

Interestingly, we observe that the rates with k = 2T correspond to nothing but those 
obtained by the thermal average. To see this point, let us introduce the averaged 7^- by 



d 3 k 
(2^)3 



d 3 k 

J2^y 



P eq (k) 



N c N D h 2 t T 

1287T 3 



(36) 



and hence (7^) 
of No 
Eq. i 



. Note that 7^(T) = (7^), which is accidentally obtained because 
2. Therefore, the assumption (ii), i.e., the thermal average of the destruction rates in 
can then be achieved by the replacement j%(k) — » (7^). It should be noted that the 



1 K 



averaged T^ q (k) in Eq. (TTTj) is given by 

K eq ) = Hl d N ) [F*F] = ^ T [F*F] , (37) 

which is exactly the (total) destruction rate of sterile neutrinos introduced by Ref. [I] and used 
in Refs. |Zl HH 023- 

In the end, by applying the assumptions (i) and (ii), our kinetic equations lead to the 
equation for Rn 

= -i[(H N ), R N ] - ^(j d N ) {F^F,R N - 1} + 2 ( 7 ^) F\A - 1)F 

-l(l d N ){F\A- l -l)F,R N }. (38) 

We should compare the equations (j551) and fl34|) with those in the previous works. Notice that 
they satisfy Eq. (135]) and hence the lepton number is conserved for each mode with k. 

In Ref. [1] the evolution of R^ and R^ has been studied by using Eq. f )38|) where only the 
first two terms are taken into account and include no effect of chemical potentials. Ref. [7] has 
solved the coupled equations of R Nj n and chemical potentials (the equations for p L i in [7] are 
nothing but those of chemical potentials), where the third term in ( 1381) and the first two terms 
in (15T]) are taken into account, however coefficients are different from the present ones. See the 
discussion below. The last terms in Eqs. f l38|) and (15T]) are original ones in the present paper. 

We shall discuss in detail the differences from Ref. [TT] where SU(2) degrees of freedom of 



lepton doublets are taken into account As for the equation (138]) there are two differences. 
One is the coefficient of the third term is "2" rather than "3" . The other one is the last term is 
missing in Ref. [TTJ. (Correspondingly, the equations for chemical potentials are also different.) 
We observe that the equations in Ref. [TTJ can be reproduced if one uses the destruction and 
production rates from the scattering process (A) or (B) (the rates for these two processes are 
the same as shown in Sec. [3]), and multiplies the rates a factor of three to get the total rates. 
This is the assumption (iii) listed before. 

As explained in Sec. [31 the processes (A) and (B) lead to the communication terms between 
Pn and pl, which induce the third term of Eq. (138]) . On the other hand, the process (C) connects 
between p N and pi and gives the fourth term of Eq. (1551) . Note again that these communication 
terms arise as the corrections to the destruction and production rates of particles, not by hand. 
As a result, we can find the explicit form of such a term with a definite coefficient. 

Having specified the differences between various sets of the kinetic equations to now, we can 
compare quantitatively the yields of the BAU. The results are summarized in Tab. [TJ where 
we take the same parameter set described in Sec. HI but both Mn = 10 GeV and 100 MeV are 



# 3 The kinetic equations in [TT] are reproduced from those in [3 [12] by replacing Tl — > Tl/Nd and fi a — > 
Nol^a- The change of Tl leads to the change of coefficient of the communication term in the equation of pN 
accordingly. The hnal value of the BAU in [TT] decreases roughly by a factor of two due to these modifications. 



Case 


I 


II 


III 


Kinetic eqs. 


Ref. HQ 


Eqs. (JMD and (1581 


Eqs. (H2D and flU} 


Momentum dep. 


Approx. (i), (ii), (iii) 


Approx. (i), (ii) 


Full 


Rates 


= (3r d ( A )) 


r<^ = /p<i(A+B+C)^ 


= r^( A+B+c )(fc) 


Y B (10 GeV) 


3.8 x lO" 10 


2.7 x 1(T 10 


2.9 x 1(T 10 


Y B (100 MeV) 


5.8 x 10" 14 


3.9 x 1Q- 14 


5.4 x 10~ 14 



Table 1: Comparison ofY B for cases I, II and III when = 10 GeV and 100 MeV. See the 
details in the text. 



considered. The cases I and II correspond to Ref. [TT] and to Eqs. ( 134)) and (|38|) . respectively. 
For comparison, we also show the results from Sec. H] as the case III. 

Let us first make the comparison between the cases I and II in which the momentum 
dependence in p N is approximately included. The momentum distributions of [pjv(fc)]22 and 
[ApAr(/c)] to t are shown in Fig. [8) Since the momentum dependence of pN is taken as (|32|) . only 
the normalization is important. As seen in the occupation numbers, the yields of N2 (and A3) 
are almost the same because the production dominantly occurs through the second term in 
Eq. (|38|) which is common for both cases. On the other hand, the asymmetries are mainly 
generated through the communication terms of the kinetic equations, and then Y B can vary 
depending on the structure of these terms. We find that the BAU in the case II is smaller than 
I by ~ 2/3, and accordingly the distributions of [Ap N (k)] to t are different by the same amount 
since Y B oc Y ANtot {= ~Y ALtot ). 

The origin of the factor 2/3 can be understood as follows: Active flavor asymmetries Y A i a 
are almost the same between two cases. This is because the three scattering processes equally 
contribute to the generation (as long as chemical potentials are small enough). On the other 
hand, only the processes (A) and (B) give a dominant contribution to Y A i tot and Y A N tot through 
the third term in Eq. (I3"8l . The contribution from the process (C) is suppressed by 0(F 2 ). 
The difference in Y B is just coming from the change of the coefficient of this term. This is the 
impact of the new communication terms found in this analysis. 

Finally, we turn to discuss the difference between the cases II and III in order to estimate 
the importance of the momentum dependence in p^. Note that the difference lies only in the 
treatment of the momentum dependence between the two cases. When we take the smaller value 
of Majorana mass as Mn = 100 MeV neutrino Yukawa coupling constants become sufficiently 
small that all the modes of interest can be away from the equilibrium one. As mentioned 
in Sec. HI the momentum dependence in the rates boost (or disturb) the production of the 
occupation numbers and also the asymmetries for the modes with lower (or higher) momenta. 
These behaviors can be seen in the distributions in Fig. [HI In this case Y B is enhanced by about 
40% due to the inclusion of the momentum dependence. 
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Figure 8: Comparison of momentum distributions of [/>#(/&)] 22 (left panel) and [Apjv(^)]tot (right 
panel) for the cases I, II and III in Tab. d The distributions are evaluated atT = Tw and are 
multiplied by (k/T) 3 . Here we take M N =10 GeV and 100 MeV. 



On the other hand, when = 10 GeV, the large Yukawa coupling constants lead to 
another effect. The low momentum modes are likely to gets in thermal equilibrium and the 
matrix of densities becomes — > p cq l. The asymmetries carried by such modes receive the 
wash-out effect. This point is clearly seen by comparing the distributions [Ap N (k)} tot of the 
cases II and III. As a result, the enhancement factor of Yb is reduced to about 10%. 



6 Conclusions 

In this paper we have investigated baryogenesis induced by the flavor oscillation between sterile 
(right-handed) neutrinos N 2 and N 3 in the i/MSM. We have presented the kinetic equations for 
the matrices of densities of sterile neutrinos pn and pfj and chemical potentials of active (left- 
handed) leptons p Va shown in Eqs. (lT2j) and (T2UI) . By using these equations the time evolution 
for the matrices of densities can be traced for each mode, and the momentum distributions can 
be found at any moments relevant for baryogenesis. As a result they allow us to estimate the 
yield of the baryon asymmetry more precisely. 
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We have evaluated the destruction and production rates of sterile and active leptons for the 
kinetic equations. Especially, we have paid a special attention to the following two respects. 
The first is to include the momentum dependence of the rate correctly, which is mandatory 
to obtain the equations of Pn,n f° r a given momentum. The second is to include the effect 
from the deviation of the initial or final state from the thermal equilibrium in the scattering 
process. As explained in Sec. [31 the second point is crucially important in the construction of 
the kinetic equations for baryogenesis in the i/MSM. We have shown that the terms connecting 
between sterile and active sectors, which are required for the successful baryogenesis, arise in a 
consistent way from the second effect to the rate. 

To be concrete, the terms connecting the matrices of densities pn and pl, which introduced 
in Ref. [7], appear as the corrections to the production rates ST P N and 5TP via the scattering 
processes (A) and (B). As shown in Sec. [5] the prefactor is different by 2/3 to the previous 
result in Ref. [TTJ. Moreover, we have found the new terms which link p^ to pi (rather than 
Pl). They are induced as the corrections 5Y d N and 5T d to the destruction rates via the process 
(C). Importantly, such a term is proportional to the product of p N ft and p Ua . Thus, our kinetic 
equations are no longer linear in Pn,Pn and p Va . It should also be stressed that both types of 
these terms are required to ensure the conservation of the total lepton number. The impacts 
of these terms on the generation of the baryon asymmetry have been described in detail. 

We have then investigated the numerical solutions of the kinetic equations. Because of 
the momentum dependence in the rates the production of sterile neutrinos is enhanced for 
the modes with k/T < 1 while diminished for the higher momentum modes. Consequently, 
the momentum distributions of p^ and p^ are significantly distorted from the equilibrium one. 
This clearly shows the importance to utilize the momentum-dependent kinetic equations. It has 
also been found that the low modes get in the thermal equilibrium and the asymmetries carried 
by them receive the wash-out effect at higher temperatures. In our choice of the parameters, 
the mode with k/T ~ 2 is the most significant for the occupation numbers as well as the 
asymmetries including the BAU. 

The comparison with the previous works has also been discussed. We have first explained 
what kind of approximations are needed to reproduce the previous kinetic equations staring from 
the present ones. It is shown that the use of the destruction rate [4] obtained by the thermal 
average is crucial to ensure the lepton number conservation when we treat the momentum 
dependence of pN approximately. We have then compare quantitatively the yields of the BAU. 
Apart from the correct treatment of the momentum dependence, the present equations gives 
the smaller BAU by a factor of ~ 2/3 compared with Ref. [11] . Within the specific choice of 
parameters, the inclusion of the momentum dependence enhances the BAU by about 40% and 
10% for M^r = 100 MeV and 10 GeV, respectively. This is because the enhancement in the 
generation of the asymmetry from the lower modes overcomes the suppression induced by the 
higher modes. On the other hand, such an amplification may be disturbed by the wash-out 
effect for large neutrino Yukawa coupling constants. 



10 



The present analysis shows that the momentum distribution of the matrices of densities 
differs from the equilibrium one and the inclusion of their momentum dependence modifies 
the yield of the BAU in non-trivial manner depending on the parameters of neutrino Yukawa 
coupling constants and masses of sterile neutrinos. Therefore, it is important to explore the full 
parameter space of the z/MSM accounting for the observed BAU by using the kinetic equations 
presented in this paper. This issue will be performed elsewhere. 
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